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It is shown that the strong volume-dependence of the axial charge of the nucleon seen in lattice 
QCD calculations can be understood quantitatively in terms of the pion-induced interactions be- 
tween neighbouring nucleons. The associated wave function renormalization leads to an increased 
suppression of the axial charge as the strength of the interaction increases, either because of a 
decrease in lattice size or in pion mass. 
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The axial form factor of the nucleon, Ga{Q 2 ), encodes 
details of the nucleon's non-perturbative structure and 
plays a key role in its properties under chiral symme- 
try p. While the axial charge, gA = Ga(0), has been 
known accurately for many years through measurements 
of neutron /3-decay, and the shape is well described by 
a dipole form with mass parameter around 1 GeV, the 
calculation of these properties within QCD still presents 
significant challenges. In particular, the mean square ra- 
dius of the axial form factor in modern lattice QCD sim- 
ulations is almost a factor two too small compared with 
experiment [2 and gA itself exhibits a remarkable depen- 
dence on the size of the space-time lattice used [3HS] ■ 

The resolution of these challenges is made especially 
urgent by the recent contributions of lattice QCD to- 
wards unravelling the origin of the spin of the pro- 
ton pi HDHT^] . For example, the isovector combination 
of the orbital angular momentum carried by the quarks, 
L u — L d , is dramatically different in several widely used 
models [I3HT5] and it can be deduced from a lattice cal- 
culation of J" — J d by subtracting the value of gA implicit 
in that same calculation. While such estimates have al- 
ready been published, one cannot assess their reliability 
without a better understanding of the systematic error 
associated with the finite lattice volume [16] . 

In this Letter we explain the origin of the dramatic 
volume dependence of gA in terms of the pion exchange 
force between neighbouring nucleons on the periodic lat- 
tice used in lattice QCD simulations. This quantita- 
tive, model independent explanation provides a frame- 
work within which to fully assess the systematic errors 
associated with the determination of the quark spin and 
orbital angular momentum as discussed above. 

One of the earliest attempts to understand why gA 
might decrease rapidly at low pion mass as the lattice 
size decreases relied on the fact that the axial current 
has a contribution proportional to V0, with <f> the pion 
field. In a finite box one may expect that the axial charge, 
given as the volume integral of the matrix element of the 
axial current would receive a contribution from a surface 
integral of the pion field, arising through Green's Theo- 



rem [17]. Promising as this idea seems, it has been proven 
that with the periodic boundary condition imposed on 
the quark fields, this surface contribution must vanish 
and hence cannot be the source of the problem [TS] . 

Next one might think of the volume dependence of the 
chiral loops which renormalize the axial charge |19j , with 
wave function renormalization reducing gA and the ver- 
tex renormalization, especially that associated with the 
axial current acting on an intermediate A baryon, tend- 
ing to increase it [2U1 [H] . These two contributions tend 
to effectively cancel each other over a wide range of bag 
radii, so that the renormalized and un-renormalized axial 
charges are very close to each other. Given that the vol- 
ume dependence of the wave function and vertex renor- 
malization is very similar for pions which do not travel 
outside the lattice volume, this is also not the source of 
the rapid variation of gA- Similar conclusions have been 
reached within other approaches [23J . 

What has not previously been recognised, however, is 
the effect associated with the fact that on a finite lattice 
each nucleon must interact with its neighbours through 
pion exchange. (Alternatively, the pion travels outside 
the lattice volume or "round the world".) In terms of 
its effect on the axial charge this may be treated in an 
analogous way to the renormalization of an individual 
nucleon. Once again the wave function renormalization, 
that is the probability to find a bare nucleon rather than 
one which has emitted a pion that will be absorbed by 
one of its neighbours, reduces gA, while the vertex renor- 
malization acts to increase it. In this case, however, 
the N — > Att transition is suppressed by an extra fac- 
tor of exp[— SL], with d — — tojv and L the lattice 
size, which is the same as the separation between near- 
est neighbours on the lattice. As a result, at all but the 
very smallest lattices (where the neighbouring nucleons 
would overlap anyway), this contribution is not effective 
at countering the effect of wave function renormaliza- 
tion. Similarly, the vertex renormalization involving Ntt 
intermediate states is small for a single nucleon (namely 
that the interchange of a z r^ with the err factors in the 
pion loop generates a factor of 1/9 compared with the 
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FIG. 1. Schematic illustration of pion emission and absorp- 
tion on a periodic lattice. 



wave function renormalization) . Hence the wave func- 
tion renormalization dominates and gA is suppressed as 
the factor exp[— m r L] increases, or m v L decreases. 

The emission of a pion which is absorbed on a neigh- 
bouring nucleon is illustrated schematically in Fig. [T] for 
a nucleon on a periodic lattice. Indeed, when the nucleon 
under consideration emits a pion, every copy also emits a 
meson and the particular one in the box under considera- 
tion must absorb a pion emitted by the other neighbour. 
Looking at the contribution to the wave function renor- 
malization coming from a single pair of neighbouring nu- 
cleons - as in Fig. [2] - one sees (in the heavy baryon ap- 
proximation which should be perfectly adequate for this 
purpose) that 
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where a, /?, a 1 and /?' represent the nucleon spin states. 

We suppose that the proton spin points up, along the 
z-axis, and consider first the interaction through the ex- 
change of a 7T° with the nearest neighbour along the y- 
axis. This is illustrated by Fig. [2] (a), where the pion 
emission necessarily flips the spin of the proton. In this 
case (T2iTu) = +1 and the contribution to the wave func- 
tion renormalization becomes 
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Defining the usual spin raising and lowering operators 

0± = \( (J x ± i<T y)' a o = a z > ( 3 ) 

we see that 

<7ly <J2y = 2 (cr 1+ (J 2 - + <Tt- °2+) , (4) 

so the spin matrix element gives a factor of 2 in this case. 
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FIG. 2. Diagrams which contribute to the wave function 
renormalization. 



Finally, evaluating the momentum integral we find 
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where K n are the modified Bessel functions (e.g., 
K (z) — !> (ir/2) 1 / 2 cxp(-z) as 2 ^ oo). 

Because of the periodicity of the lattice one can also 
exchange a tt + with the neighbouring nucleon along the 
y-axis. In this case the isospin factor is two, rather than 
one, and hence 8Z [ 2 h) = 28Z ( 2 a} . For the neighbouring 
nucleon along the z-axis the spin factor is unity. That is, 
the spin is not flipped. Thus SZ 2 and 5Z 2 are each half 
as large as SZ^ and 8Z^ , respectively. This asymmetry 
between lattice sites perpendicular and parallel to the 
nucleon spin is a natural consequence of the tensor force 
generated by the exchange of a pseudoscalar meson. 

The expressions above give the contribution to the 
wave function renormalization associated with pion ex- 
change between a single pair of nucleons. To ensure that 
we represent what goes on on the lattice we need to 
include all neighbours which make significant contribu- 
tions. There are four nearest neighbours in the x— y plane 
and two along the z direction. Thus the total contribu- 
tion to the wave function renormalization from nearest 
neighbours, SZ 2 n , is 
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As explained earlier, for gA the vertex renormaliza- 
tion is not very effective in countering the wave func- 
tion renormalization when the A excitation is suppressed. 
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However, for a lattice which is not spherically symmet- 
ric the suppression is not simply 1/9 but differs from 
nearest (nl) to next-to-nearest (n2) to next-to-next-to- 
nearest neighbours (n3). It is a trivial calculation to show 
that for nl the vertex renormalization is one fifth of SZ™. 
Thus the total relative correction to gA (he. .9a (1 — <5<m)) 
from nearest neighbours is: 
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FIG. 3. Illustration of the convergence of the correction to 
gA calculated on a lattice 3.5 fm on a side as we include the 
nearest neighbour (nl) up to next-to-next-to-next-to-nearest 
neighbour (n4), for three values of m^L. 



Turning to the effect of the next-to-nearest neighbours, 
we note that there are more of them (12 instead of 6) but 
that their distance from the nucleon under study is y/2L. 
Thus, provided m^L is relatively large, we expect that 
the correction to gA should converge relatively quickly. 
In the present work we consider the corrections out to 
and including next-to-next-to-nearest neighbours, or a 
distance \^3L. In these cases, the factor one fifth for nl 
becomes one sixth for n2 and one ninth for n3. In sum- 
mary, the total corrections to gA from next-to-nearest 
and next-to-next-to-nearest neighbours are 



dgf = -24T(ml,V2m^L) 
o 



Sgf = -12F(ml,V3m v L) 



The degree of convergence for three relatively light quark 
masses is illustrated in Fig.|3]for a 3.5 fm lattice. However 
we point out that the convergence of the correction is 
solely dependent on to^L, and the figure displays the 
convergence down to m„L ~ 2.5. Here we also show the 
contribution from the next (2L) term which justifies our 
truncation at a distance of V3L in the remainder of this 
letter. 

The total reduction in gA to this order is illustrated in 
Fig. [3] — note that for each volume, the curves are plotted 
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FIG. 4. The total reduction in gA from wave function and 
vertex renormalization is illustrated as a function of pion mass 
(in GeV) for several lattice sizes. 



only down to a minimum reliable pion mass estimated 
by m^L — 2.5. This figure makes it clear that these 
corrections are substantial. Modern lattice simulations 
of gA are now available at masses as low as 0.3 GeV, but 
even there the correction is almost 20% on a 2 fm box and 
6% at 2.5 fm. Even more significant, we see that as the 
pion mass approaches the physical value the correction 
is as large as 10% on a 4 fm lattice. This makes very 
clear the challenge of a brute force determination of this 
quantity at the physical quark mass. 

We show the effect of this correction applied to a sim- 
ple fit linear in m%, g A (l — <5#a) + Brn^, with two fit 
parameters (g A ,B) applied to recent lattice data from 
the QCDSF collaboration in Fig. [5] and the RBC collab- 
oration in Fig. [6] Note that we have made no attempt 
here to include any non-analytic effects in the pion mass 
dependence of gA, although in this case, these effects are 
anticipated to be small [H Q]5] . The very dramatic reduc- 
tions in gA at small lattice sizes and low pion mass are 
very clearly seen there, and we find remarkable agreement 
between the present calculation and those data. 

Earlier work has suggested that similarly large finite- 
volume effects can be described by the discretisation of 
standard one- loop corrections [3] . The enhancement seen 
there can by understood noting that both the input g&N 
coupling and the fitted <?aa coupling are about 1.5 times 
larger than those predicted by SU(6) relations. In the 
second panel of Fig. [5j the standard contributions using 
SU(6) estimates for the couplings are shown for the two 
relevant volumes. We see that the effects are insignif- 
icant compared to the exchange effects investigated in 
the present work. 

The calculations reported here show clearly that a sig- 
nificant part of the hitherto unexplained quenching of 
gA on small lattice volumes originates with the pion ex- 
change force between neighbouring nucleons on a periodic 
lattice. In retrospect, this explanation seems to be very 
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FIG. 5. The axial charge gA, with the finite volume correction 
associated with pion exchange between neighbouring nucleons 
included, is plotted for several different values of the lattice 
size, L (in fm). Also plotted are lattice calculations of qa 
from QCDSF [U[9] with three different /3's (lattice spacings), 
/3 = 5.25 (a = 0.076 fm) /3 = 5.29 (a = 0.072 fm) /3 = 5.40 (a = 
0.060 fm), on multiple volumes. For comparison, in the middle 
figure we also show (dot-dash line) the predicted finite size 
effects from g], using the SU(6) A couplings |22j . 
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FIG. 6. Curves as described in Fig. [51 plotted with lattice 
calculations of gA by Yamazaki et al. [3]. 



natural, although the key role of wave function renormal- 
ization was somewhat unexpected. 

With this fascinating issue resolved, one may ask what 
other problems being addressed by lattice QCD calcula- 
tions may be similarly affected. The study of the fraction 
of the nucleon spin carried by its quarks is certainly a top- 
ical example, as is the orbital angular momentum of the 
quarks and gluons and hadron magnetic moments. The 
tensor nature of the pion exchange force suggests that 
special care may be required in those cases where spin 



related deformation is important [24 , for example in the 
study of the E2/M1 ratio for the A. 
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